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1 Introduction 

Both the most challenging and the most important question one needs to address when 
trying to solve numerically a Coefficient Inverse Problem (CIP) is: How to obtain a rigor¬ 
ously guaranteed good approximation for the exact solution without any advanced knowl¬ 
edge of a small neighborhood of this solution? We call a numerical method addressing this 
question globally convergent. The reason of the importance of this question is that con¬ 
ventional least squares cost functionals for CIPs are non-convex. Hence, they have many 
local minima and ravines. This makes numerical procedures, which use these functionals, 
unreliable. Indeed, to find a proper minimum, one should start the minimization process 
in a small neighborhood of the exact solution, i.e. one should start from a good first guess 
for the solution. However, if such a guess is known a priori, then why the exact solution 
is also not known a priori? In addition, such a neighborhood is very rarely available in 
practice. Therefore, it is important to develop globally convergent numerical methods for 
CIPs. 

Currently there exist two types of such methods. Both for the case of the non overde¬ 
termined data, i.e. for the case of a single measurement event. The first type is based 
on constructions of globally strictly convex Tikhonov-like functionals. Carleman Weight 
Functions (CWFs) are the key to the global convexity. This method was initiated in 
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works of Klibanov in 1997 mm3; also see Klibanov and Timonov |13j . Recently there 
is a renewed interest in Beilina and Klibanov [7] and Klibanov and Thanh [T5] with some 
numerical studies in B3j. A different, although a similar approach, was carried out by 
Baudouin, De Buhan and Ervedoza [2j. The global convexity is understood as follows: 
Given a convex set G of an arbitrary diameter d in a certain Sobolev space, one can choose 
the parameter Ao = Xo(d) » 1 of the CWF such that for all A > Ao that functional is 
strictly convex on K. Assume now that there exists a minimizer of that functional on the 
set K. Then the strict convexity guarantees convergence of the gradient method to this 
minimizer starting from any point of the set K mug. This is the global convergence as in 
the above definition. Recall that in the conventional case the gradient method converges 
to a minimizer only if starting in a small neighborhood of that minimizer. 

The second type of globally convergent numerical methods for CIPs is the method, 
which was initiated in the paper of Beilina and Klibanov [3] and was discussed since then 
in a number of follow up publications of these authors with coauthors. Results obtained 
before 2012 were summarized in the book [Hj. In particular, this method was completely 
verified on experimental data, see, e.g., Chapter 5 in [4j, as well as mm- We also refer 
to a recent paper of Chow and Zou [9] for this method. 

The method of this paper falls into category of the first type of globally convergent 
numerical methods. Using a CWF for the parabolic operator, we construct a globally 
strictly convex cost functional for a CIP for a general parabolic equation of the second 
order. Unlike this, Note that CIPs for hyperbolic PDEs were considered in [2J j7[ [IT, 
ng. Although a CIP for a parabolic PDE was considered in |12j . the main difference of 
that work with the current paper is that in [12] a certain series was truncated, whereas 
truncation does not take place here. 

In section 2 we state our inverse problem. The globally strictly convex cost functional 
for it is constructed in section 3. In section 4 we prove the main result of this paper, 
which is Theorem 1. 


2 Statement of the problem 


Let 12 C K 3 be a bounded domain with a piecewise smooth boundary dfl. Let T > 0 be 
an arbitrary number. Denote = Ox (— T, T ). Let T C dVt be a part of the boundary 
dVt and let T eh 2 . Denote = Y x (— T, T ) . Consider the elliptic operator of the second 
order in D, 

3 3 

Lu=^ 2 a ij ( X ) U XiXj + Ux i + c ( x ) u , x e (2.1) 

i,j =1 i,j =1 

Here 

= aji , G C 1 (O) ; b j: c G C (JT) , (2.2) 
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Pi lei 2 < E (*) tej ^ ^ l^l 2 > V * G n, G ID, (2.3) 

i,j= 1 

where p 1; p 2 — const. > 0,p 1 < p 2 - Let the function u G C 4,2 (Or) satisfies the following 
conditions 

u t = Lu in fl T , (2.4) 

u{x,0) = f(x). (2.5) 
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We assume here that equation (12.41) is valid not only for t > 0, but for t < 0 as well. 
This is because the Bukhgeim-Klibanov method [5] does not work for the case when a 
parabolic equation is valid only for t > 0 EEJ Ej • Our interest is in the inverse problem 
which we now formulate. 


Coefficient Inverse Problem (CIP). Assume that the coefficient c(x) in H2.1\) is 
unbiown for x G D. On the other hand, assume that the following functions g\ (x, t ), g 2 (x, t ) 
are known, 


U |r T = 9i (x, t ), d n u |r T = g 2 (x, t ). 

(2.6) 

Determine c (, x) for x G O. 


Assume that in (12.51) 


/ ((c) >2 b = const. > 0 in D, d . 

(2.7) 


Then uniqueness of this CIP follows immediately from Theorem 1.10.7 of [4j as well as 
from Theorem 3.3.2 of [T3] and Theorem 3.4 of [Ej. Note that a particular case of the 
data (]2.6I) is the case of backscattering data. The principal parts of the operators L and 
d t — L we denote as L 0 , P 0 , 
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Lqu ^ ^ aij (x) u XiX j, PqU Ui Lqu. (2-8) 

m=i 


3 The cost functional 


To construct a globally convergent numerical method for our CIP, we construct in this 
section a globally strictly convex cost functional. Below x = (xi,x) , where x = (x 2 ,xf) . 
Without any loss of generality we can assume that T = j;r : x\ = p (x ), p G C 2 ^|T| < Vd)} 
for a number d G (0,1). Therefore, changing variables as (xi,x) 44 ( x[,x ), where 
x[ := x\ —p(x), and using the same notations as before, for brevity, we conclude that we 
can assume that 

T = [x : x — (0, x) , |x| < VdX . (3.1) 


Denote 


Q =< x : xi > 0, \x\ < Vd 


(3.2) 


and assume below that f Y C D. Then T C dfl d . Below we determine the unknown 
coefficient c ( x ) only in a subdomain of the domain Q d . 

We now formulate the Car Ionian estimate for the operator Pq = dt — L$. Let A > 1 
and v > 1 be two large parameters, which we define later. Consider an arbitrary number 
a G (0, d). Consider functions ( x , t), ( x , t) defined as 
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if (x, t)=x 1 + \x\ 2 + — + a, ip x (x, t ) = exp . 

Consider the following sets 


G a , d = <! (x, t) : X! > 0, xi + |x| + — + a < d 


Ga,d ~ G a ,d n {t — 0} , 
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(3.3) 


L,d = {(x,t):x 1 >0,x 1 + \x\ 2 + J— + a = d}, 
r a,d,T = <J (x, t) : x 1 = 0, \x\ 2 + ^<d-fl[cr T , 


dG a . d — U r 0)d)T , 


Ga,d-e = <j ( X , t) : X 1 > 0, £i + |x| 2 + ^ + a < d - £ \ , 


where e G (0, d — a) is a sufficiently small number. Clearly G a ^- £ C Since d G (0,1), 
then G a) d C fir, G° ad C Xl d Also, f a is the level surface of both functions f>, <^ A . Note that 
f ad is the level surface of the function ip x , 


ffiin (fix = vl k a , d = exp [2A d v ] . (3.4) 

Ga,d 


Lemma 1 follows immediately from Lemma 3 of §1 of chapter 4 of the book of Lavrentiev, 
Romanov and Shishatskii [ IB] . 

Lemma 1 (Carleman estimate). There exist sufficiently large numbers u 0l Aq > 1, 


uq uq ^a, d, /i 1; p. 2 , max lkj|| c i^ , ? Ao Ao ^n, d : P 21 m.cpc ,T 

depending only on listed parameters and a sufficiently large absolute constant Ao > 1 such 
that for all v > v 0 , A > A 0 and for all functions u G C 2,1 (G atd ) the following pointwise 
Carleman estimate is valid for all (x, t ) G G a ^ 

(P 0 u) 2 (fl > CX |Vu| 2 + CX 3 u 2 (p\ + div U + V t , 

\u\,\v\ < CX 3 [(\7u) 2 + u 2 t +u 2 ]vl 


where the constant C 
Denote 


L c u 


C 


Tl, HlQX.ij 


> 0 depends only on listed parameters. 
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Lu — C (x) U = ^2 a ij ( X ) u XiXj + ^2 ( x ) u *i ■ 

i,j=l i,j =1 


Since the function u G C 4,2 (fir) , then (12.711 implies that u(x,t ) > b > 0 in G a)d for 
sufficiently small T. Hence, we can consider the function v (x, t ) = In u (x, t ). Substituting 
u = e v in (12.41) and (12.51) , we obtain 


v t = L c v + ^2 a i,j ( x ) v Xi v Xj + c (x) in G a , rf , (3.5) 

i,j =1 


v (x, 0) = In / (x ). (3.6) 

Let w (x,t) = v t (x,t ). Differentiate (13.5|) with respect to t and use (13.61) . We obtain the 
following nonlinear integral differential equation with respect to the function w in the 
domain in G a , d 

w t = L c w+ 
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+ Y1 a 0J ( x ) w *i ( ( ln f)xi + / w *3 r ) dr 
i,j =1 


(3.7) 


+ Y1 ai ’i ^ Wx i ( ln f)xi + / w *i (®» r ) dT 

ij=i \ o 

In addition, conditions (I2.6[) imply that 

w |r„, Ar = 9i (x,t),d n w Ir.„= 92 OM), 


(3.8) 


where 


9i (x, t) = d t ln^i (x, t) , 52 (x, t) = — ~ 

9 1 Pi 


Since functions 51,52 are the data for the inverse problem, then they naturally contain 
noise. Even though the differentiation of a noisy function is an ill-posed problem, it can 
be handled by a number of well known regularization methods, see, e.g. Aristov [tj. 

Thus, we have obtained the nonlinear integral differential equation (13.711 with the 
lateral Cauchy data (13.811 . If we find the solution of this problem, then backwards calcu¬ 
lations will deliver us the target coefficient c (x) for x G G° ad . Hence, we focus below on 

the solution of the problem (13.711 , (13.811 . Denote Lw the right hand side of (13.711 . plus the 
term —w t , 

Lw = —w t + L c w+ 

+ a D' (A) w *i ( ln f)x d + / w *j r ) dr I ( 3 - 9 ) 

0=1 \ n 


+ <Hj (®) w x . | (In f) x . + j w Xi ( x , r) dr 

0=1 


Our weighted Tikhonov-like cost functional is 


J\,a (w) = exp ^— 3Xd l/ ) / yLwj ip\dxdt + a\\w\\ H4 / G ^ 


(3.10) 




where a G (0,1) is the regularization parameter. We use the multiplier exp (—3A d~ u ) to 
ensure that we can indeed choose a G (0,1), see Theorem 1. We use the H 4 ( G a; d ) — norm 
here since we need in our proof w G C 1 ( G a) d ) fl H 2 (G a ,d) , and the embedding theorem 
guarantees that 


H‘ (G a4 ) c C 1 (G m ) , ||«|| ci(3b d) 


<C\ 


u 


H 4 (G aid ) 


, Vm € H“ (G m ), 


(3.11) 


where C\ = C\ ( G a> d ) > 0 is a generic constant depending only on the domain G a ,d- Thus, 
we consider below the following problem. 

Minimization Problem. Minimize the functional J\ tDl (w) in US. 10 1) . subject to the 
lateral Cauchy data \3.8\) . 

Let R > 0 be an arbitrary number. Consider the set B (R) C H 4 (G ai d ), 


B(R) = 
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e H 4 (G a , d ) : \M Hi , Ga d ^ < R,w satisfies boundary conditions eg}. (3.i2) 

Introduce the space Hq (G ad ) as 

Ho (G a , d ) = {u G H 4 (G a:d ) : u \r atd}T = d xi u |r a , d , T = 0} . 

Theorem 1. For all numbers A, v, a > 0 and for all functions w G B ( R ) there exists 
the Frechet derivative J'\ a (w) G (G a ,d) of the functional J\ t0l at the point w. Let 
v = Vo be the sufficiently large number of Lemma 1 and b > 0 be the number in \2. 7p . 
There exists a sufficiently large number Ai, 


Ai = Ai {a, d, ix 2 ,max\\ a i,j\\ c i(p^ ,T,b,R, ||V/|| c ^o^J > A 0 > 1 (3.13) 

depending on listed parameters such that if the regularization parameter 

a G (exp (—A/ (2cT)), 1), then the functional J\,a(u>) is strictly convex on the set 
B (R) for all A > Ai. More precisely, 

j\,a (w 2 ) - j\ a (wi) - J\.a ( w l) i w 2 ~ Wl) > 

C\ exp (2Ag) J [(Vw 2 - V%) 2 + (w 2 - Wi) 2 ] dxdt + ^ \\w 2 - Wi\\ 2 H ^ Ga d ) , (3-14) 

Ga+e,d 

Wwi, w 2 G B ( R ), VA > Ai. 

where q = (d — e)~ u ° [1 — 3 (d — e) u ° / (2^)^°] > 0 and the constant Ci > 0 depends on the 
same parameters as those in \3.13\) . 

Note that we can require a G (exp (—A/ (2d v )),l), since exp (—A/ (2 d u )) « 1 for 
sufficiently large A. This theorem is the main result of our paper. Theorem 1 enables 
one to prove the convergence of the gradient method, which can start at any point of the 
set B (. R ). Since there are no restrictions on the diameter 2 R of this set, then this is the 
global convergence as defined in Introduction. The stability with respect to the noise in 
the data 03.81) can also be established using this theorem. These two latter results can 
be derived from Theorem 1 in the same manner as similar results are derived from global 
strict convexity theorems in izma. Hence, we do not describe these results here. Below 
we focus on the proof of Theorem 1 and assume that its conditions are satisfied. Below 
C\ > 0 denotes different constants depending on the same parameters as ones listed in 

(HDD- 


4 Proof of Theorem 1 

Let W\ , W ‘2 G B (R) be two arbitrary functions. Denote h — w 2 — W\. Then 

h 6 H* (G m ) , < 2R. (4.1) 

Let A = (^L (wi + /i) j — ^L (uq) j . First, we single out the linear part of this expression 
with respect to h. Using (12.ip and (13.91) . we obtain 

L (wi + h) — L (-uq) + L c (h ) — h t 
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+ Y a i,j ( x ) K (In /) + / W lxj (x, r) dr + 


+ Y ai ’i h *i ( ln f)xi + / w ^i r ) drT + 


Hence, 


+ Y ai ’i ( x ) Wl ^ / h xj ( x ’ r ) / h x t r ) • 


A = 2L (wi) L(h) - ht+Y a tj ( x ) h *i ( ln /) Xj + / (®, r) dr 


+2 L (wi) ^ a, ;j (x) h Xj ( (In f) x . + / Wn, (x, r) dr 


+2L (wi) y; a i:j (x) w lxi h Xj (x,r) dr+ w lxj h Xi (x,T) dr + S 2 (h,w i), 


where 


S' (h, wi) = L(h) — h t 

3 / y 

+ X] ai ’i ( x ) h *i ( ln /)*,- + / (a, r) dr 


+ aij (x) h Xj (ln f) x . + / (x, r) dr 


+ Y ai ’i ( X ) Wl Si / h Xj r ) drr + / h Xi (x, T ) d,T 


The expression D (uq, h) = A — S 2 (h, w\) is linear with respect to h. Hence, consider the 
functional 

Qa, q (h) = exp (— 3Xd~ u ) / D (wi, h ) c p\dxdt + 2a [/i, wi], 

G a , d 

where [, ] is the scalar product in H A (G a> d) ■ This is a linear bounded functional acting 
from Hq (G a ^) into R. Hence, by Riesz theorem there exists unique element U\ }01 (uq) G 
Hq ( G at d ) such that Q\ tCt {h) = [U Xa (w\) , h]. Furthermore, the norm \ equals 

to the norm of the functional Q\, a - Hence, we have proven the existence of the Frechet 
derivative J Xa (uq) = U Xa (wi) G Hq (G a ^) of the functional J Xa and 


J X a (wi) ( h ) = exp (—3A d u ) / B (uq, h ) <p\dxdt + 2a [h, uq], \/h G Hq (G a ,d) 
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Hence, 


j\,a (Wl + h) - Jx,a (w i) - J' Xa {w i) (h) 


(4.3) 


= exp (—3Ad J S 2 ( h , uq) (p 2 x dxdt + a 

G a , d 


H*(G a}d ) ’ 


where S (h,uq) is given in (14.2[) . 

We now focus on the estimate from the below of the integral in (14.3jl . Because of 
Lemma 1, we single out the term with (h t — L(h)) 2 . Using (12.21) . (12.71) . (13.9[) and the 
Cauchy-Schwarz inequality, we obtain 


S 2 (h,w 1 )>-(h t -L(h)) 2 -C 1 (Vh) 2 -C 1 ( / | Vh (x, t)\ 2 dr 


> l -(P Q (h)) 2 -Ci{Vhf-Ci | / \Vh(x,r)\dr ) , 


(4.4) 


where the operator P 0 — d t — L 0 was defined in (12.81) . It follows from Lemma 1.10.3 of [4j 
that 

t 


| V/i (x, t)\ dr | ip\dxdt <f (Vh) 2 <p\dxdt. 


(4.5) 


G a , d \ 0 


G„. 


Hence, using Lemma 1, (13.31) . (13.4]) . (14.11) . (14. 4 p and (14.5p . we obtain for sufficiently large 

A > Ai 


exp (—3Ad j S 2 (h, uq) ip\dxdt > Ci\ exp (—3Ad ") J [(V/r) 2 + /i 2 ] Lp\dxdt 

G a ,d G a ,d 


—C\ exp (—3A d I/ ) J (Vh) < p\dxdt — CAexp (—Ad ") J [h 2 + (Vh) + h 2 ] da 

G a ,d £>a,d 

> UiAexp (—3A d~ u ) j [(Vh) 2 + /r 2 ] <p 2 x dxdt — C\ exp (—A d~ u ) J [h 2 + (Vh) 2 + h 2 ~\ da 

Ga,d £ a,d 

> C\ exp (2Ag) J [(V/r) 2 + h 2 ] dxdt — C\ exp (—A d~ u ) j [h 2 + (V/r) 2 + h 2 ] da. 

G a ,d—e £a,d 

Thus, we have established that 


exp (—3Ad ") / S 2 (h, wi) <p\dxdt 
G a , d 


(4.6) 


> Ci exp (2A q) / [(V/r) 2 + h 2 ] dxdt — Ci exp (—Ad / [h 2 + (V/r) 2 + h 2 ] da. 


G „ d _ 


^ <a,d 












Next, since a G (exp (—A/ (2d?)) , 1), then 


-C, exp (—Ad-) J [fcf + (Vfe) 2 + ft 2 ] dcr >-| ||h£ l(o>i) . («■) 

Ca,d 

Combining (14. 3 J) with (I4.6j) and (14.71) . we obtain the target estimate (I3.14p . □ 
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